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ABSTRACT 

HH 900 is a peculiar protostcllar outflow emerging from a small, tadpole-shaped glob¬ 
ule in the Carina nebula. Previous Ha imaging with HST /ACS showed an ionized 
outflow with a wide opening angle that is distinct from the highly collimated struc¬ 
tures typically seen in protostellar jets. We present new narrowband near-IR [Fe n] 
images taken with the Wide Field Camera 3 on the Hubble Space Telescope that reveal 
a remarkably different structure than Ha. In contrast to the unusual broad Ha outflow, 
the [Fe n] emission traces a symmetric, collimated bipolar jet with the morphology 
and kinematics that are more typical of protostellar jets. In addition, new Gemini 
adaptive optics images reveal near-IR H2 emission coincident with the Ha emission, 
but not the [Fe 11]. Spectra of these three components trace three separate and distinct 
velocity components: ( 1 ) H2 from the slow, entrained molecular gas, ( 2 ) Ha from the 
ionized skin of the accelerating outflow sheath, and ( 3 ) [Fe 11] from the fast, dense, and 
collimated protostellar jet itself. Together, these data require a driving source inside 
the dark globule that remains undetected behind a large column density of material. 
In contrast, Ha and H2 emission trace the broad outflow of material entrained by 
the jet, which is irradiated outside the globule. As it get dissociated and ionized, it 
remains visible for only a short time after it is dragged into the H 11 region. 
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1 INTRODUCTION 

Protostellar outflows are a beacon of star formation, signal¬ 
ing active disc accretion even in deeply embedded regions 
where the disc, and sometimes protostars themselves, can¬ 
not be detected directly. We discuss one such example in 
this paper. The detailed physics o f .jet launch and col lima- 
tion is not yet understood (see, e.g. [Ferreira et al]|2006i '). but 
disc accretion ultimately must fuel the jet. Observations of 
relatively unobscured young stars measure this jet directly, 
but trace an epoch long after the most active accretion. CO 
observations trace outflows from young sources and those 
in more embedded regions. The observed emission may be 
dominated by am bient molecul ar material entrained by the 
jet (e.g. Arcc & Sargent 2005) and/or from a slower disc 
wind (e.g. Klaassen et al.1 20131 '). 

Spatially resolved outflows can be a powerful tool to 
study star formation at higher stellar mass where evi¬ 
dence for discs remains elusive. There have been only a 
few direct detections of disks around intermediate-mass 
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protostars (e.g. Kraus etah l20ld : IPreibisch et al.l l201ll : 
ICarrasco-Gonzhlez et al.l l2012l ~). Structure in the form of 
clumps and knots along the jet axis points to the vari¬ 
able nature of accretion and outflow, and provides one of 
the few ways to infer a given protostar’s accretion history. 
The longest jets extend parsecs on the sky, sampling a 
signif i cant fraction of th e accretion age (e.g . iMartf et al.l 
1 19931 : iDevine et al.l Il997l : ISmith et alJ l2004af ). Many out¬ 
flows from moderate- to high-mass young sta rs appear to be 
scaled-up versions of low-mass sy stems (e.g. iGuzman et al.l 
l201ll : iReiter fe Smitbl 120131 . 20141'), although t his is not al¬ 
ways the case (see, e.g. Shepherd et al.l [20031 ). Other out¬ 
flows, such as the Becklin-Neugebauer / Kleinmann-Low 
(BN/KL) outflow in Orion OMC1, show evidence for vio¬ 
lent processes related to the dynamical interaction and ejec¬ 
tion of high-velocity stars that creates an explosive morphol¬ 
ogy (iBallv et al.ll201ll : iGoddi et al .11201 il l . However, similar 
outflow behavior from a large number of protostars over a 
wide mass range provides compelling evidence that mas¬ 
sive star formation can be understood as a scaled-up ver¬ 
sion of low-mass star formation. It is then essential to study 
intermediate-mass stars (~ 2 — 8 Mq) where any changes in 
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the dominant physics of formation between low- and high- 
masses are expected to occur. 

Outflows from young, intermediate-mass protostars are 
typically observed at millimeter wavelengths where emission 
from entrained molecules penetrates the high column den- 
sities that characteriz e massive star- f ormin g regions (e.g. 
Thjcahashi et al.l 120071 : iBeuther et al.l 120081 : iBeltran et al] 


. However, in an H II region, much of the obscuring gas 
and dust may have been cleared. In such environments, ul¬ 
traviolet (UV) radiation from nearby massive stars will illu¬ 
minate protostellar jets propagating into the H II region cav¬ 
ity. Unlike outflows in quiescent regions where emission at 
visual wavelengths traces only shock-heated material, exter¬ 
nal irradiation in H II regions lights up the body of the entire 
jet, revealing outflow material that would otherwise remain 
invisible. This allows the physical properties of irradiated 
jets to be measured using the diagnostics of photoioniz ed gas 
(e.g. iBallv fc Reipurthll200ll : lYiisef-Zadeh et, al .Il2()05h with 
the high angular resolution observations that are available 
at shorter wavelengths. 

HH 900 is on e of the many HH jets discovered by 
ISmith et all (I2OI0I ') as part of their Ha? survey of the Ca¬ 
rina nebula using the Advanced Camera for Surveys (ACS) 
onboard the Hubble Space Telescope (HST). HH 900 is an un¬ 
usual bipolar outflow emerging from a small tadpole-shaped 
globule located ~ 3 pc (in projection) away from r? Carinae 
see F igure CD. In earlier ground-based images. [Smith et al.1 
200311 identified it as a candidate proplyd, although nei¬ 
ther of the two ‘tails’ emerging from the globule point away 
from 77 Car. Higher resolution HST images show that these 
two ‘tails’ appear to be a wide bipolar outflow, and reveal 
a number of other peculiarities. Along the western limb of 
the broad bipolar outflow lies a strong Ha filament and a 
point source nearly at its center, raising the possibili ty th at 
the filament is a microjet driven by the star llSmith et al.l 
12010 j. Without kinematic information, it remains unclear 
whether the point source actually drives the putative Ha 
microjet, and if so, whether it is physically related to the 
larger bipolar outflow. Estimated mass-loss rates (derived 
using the Ha emission measure) from the wide inner jet in 
the eastern limb of the outflow and the microjet along the 
western limb are the highest of all the outflows found by 
[Smith et al.l ( 20ld). 


Shinn et al.l (|2013i f present ground-based observations 


of [Fe 11 ] 1.64 pm emission in the Carina nebula. Based on 
the morphology of the [Fe 11 ] emission from HH 900, they 
propose that the HH 900 driving source is one of the young 
stellar objects (YS Os) mo deled in the Pan-Carina YSO Cat¬ 
alog fPCYC. IPovich et al.ll201lh . A second Spiteer-identihed 
YSO lies along the western limb of the flow, although with 
angular resolution > 1.5”, it i s unclear how th i s sour ce re¬ 
lates to the putative microjet. lOhlendorf et alj (12012)4 con¬ 
clude that Herschel emission near HH 900 is from the ex¬ 
ternally heated globule and is unlikely to come from young 
stellar objects in the region. 

I 11 this paper, we present new optical and IR obser¬ 
vations of HH 900 that allow us to investigate its mor¬ 
phology and kinematics in detail. New narrowband IR im¬ 
ages obtained with HST probe the [Fe 11 ] 1.26 pm and 
1.64 pm lines that are often assumed to be shock-excited 
in protostellar jets, although this is not necessarily the 
case in regions with significant FUV radiation. Indeed, 


iReiter fc Smithl (l2013l i showed that these two [Fe 11 ] lines 
probe dense, low-ionization jet material not traced by Ha. 
Near-IR [Fe 11 ] lines can also penetrate the extinction inside 
the dusty birthplaces of these jets, allowing us to connect 
the larger Ha outfl ows to the gpfteer-i dentified protostars 
that drive them l|R.eiter fc Smith|l2013l f. Both [Fe 11 ] 1.26 
pm and 1.64 pm originate from the same a 4 D level, so their 
flux ratio is insensitive to excitation conditions. Variations 
in the flux ratio along the jet therefore provide one way to 
measure variations in the reddening of the immediate jet en¬ 
vironment. This may be particularly interesting in HH 900 
where extended H 2 emission (Hartigan et al. 2015, in prepa¬ 
ration) and dark filaments in Ha images suggest the presence 
of molecules and dust in the broad bipolar flow. 

We also present new, second-epoch Ha images from 
HST that allow us to measure jet motions in the plane of 
the sky. The outward motion of jet knots provides a direct 
kinematic identification of candidate jet driving sources. To¬ 
gether with ground-based optical and IR spectra, our Ha 
proper motions and [Fe 11 ] images from HST offer a new 
view of HH 900 that simultaneously help to unravel but also 
deepen the mystery of this unusual outflow. 


2 OBSERVATIONS 

2.1 High-resolution images 

[Fe 11 ]: Table |T] lists the details of the images and spec¬ 
troscopy presented in this paper. Near-IR, narrowband 
[Fe 11 ] images and corresponding off-line narrowband con¬ 
tinuum images were obtained with WFC3-IR on board the 
Hubble Space Telescope (HST) under program GO-13391 
(PI: N. Smith) in Cycle 21. For HH 900, we duplicated the 
observing strategy used previously to obtain [Fe 11 ] images 
of four other jets in t he Carina nebula, HH 666, HH 901, 
HH 902, and HH 1066 (IReiter fc Smithll2013l 'l. We employed 
a box-dither pattern to avoid dead-pixel artifacts and to 
provide modest resolution enhancement. Using the same in¬ 
tegration time for all filters provided a similar signal-to-noise 
ratio for the [Fe n] and adjacent continuum images. 

H 2 : We obtained near-IR., narrowband H 2 2.12 pm 
images and complementary off-line narrowband K(short) 
continuum images (central wavelength of 2.093 pm) dur¬ 
ing early science wi th the Gemini South Adaptive Op- 
tics Imager (GSAOI. iMcGregor et al1l2004l : ICarrasco et al.l 
l2012l l. GSAOI is a near-IR imager, used in combination with 
Gemini Multi-Con j ugate Adaptive Op tics System (GeMS, 
iRigaut et al.l 120141 : iNeichel et al.l I2014T ) with natural guide 
stars to provide near-diffraction limited images over a field 
of view of 85' x 85' in the 0.9-2.5 pm range. HH 900 was 
observed 22 March 2013 in queue observing mode with a 
9-point dither pattern and individual integrations of 120 s. 
We resampled the aligned, flat-fielded images to the same 
pixel scale as the WFC3-IR images. Simultaneous K(short) 
continuum images allow us to subtract bright IR continuum 
emission from the edge of the globule and therefore search 
for extended H 2 emission from the jet. 

Ha: Ha images of HH 900 were ob tained with 
HST/ ACS for the first time on 18 July 2005 ([Smith et al.l 
l2010h . We obtained a second epoch on 04 August 2014 us¬ 
ing the the same instrument and filter, duplicating the first- 
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Table 1. Observations 


Instrument 

Filter / Position 

Date 

Int. time 

Comment 

ACS /HST 

F658N 

2005 Jul 18 

1000s 

Ha + [N II] 

ACS /HST 

F658N 

2014 Aug 04 

1000s 

Ha + [N II] 

GSAOI/ Gemini 

K s N cont. 

2013 Mar 23 

1080s 


GSAOI/ Gemini 

h 2 

2013 Mar 23 

1080s 


WFC3-IR/RST 

F126N 

2013 Dec 28 

2397 s 

[Fe ii] A12567 

WFC3-IR/RST 

F128N 

2013 Dec 28 

2397 s 

continuum 

WFC3-IR/HST 

F164N 

2013 Dec 28 

2397 s 

[Fe n] A16435 

WFC3-IR/-H1ST 

F170N 

2013 Dec 28 

2397 s 

continuum 

EMMI/JVTT 

Ha, [S n] 

2003 Mar 09 

900 s 

P.A. = 63° 

FIRE/Magellan 

West 

2014 Jan 15 

600 s 

P.A. = 77° 

FIRE/Magellan 

East 

2014 Jan 16 

600 s 

P.A. = 60° 


epoch observational setup in order to measure proper mo¬ 
tions (program GO-13390, PI: N. Smith). Together, this pro¬ 
vides a ~ 9 yr time baseline between observations, allowing 
us to measure the motion of faint jet features to ~ 25 km 
s -1 . Using the same orientation and coordinates for both 
epochs minimizes position-dependent errors when determin¬ 
ing proper motions. 

We follo w a method similar to that of I And erson et ahl 
([2008a| 3), lAnderson fe van der Mare! ( 20ld 'l and 
ISohn et al.l (|20l2h to align and then measure precise 
proper motions in the ACS images. This method is based 
on PSF photometry of the bias-subtracted, flat-fielded, 
and CTE-corrected flc images produced by the HST 
pipeline. Unlike the drizzled drc images, which have been 
calibrated, flat-fielded, CTE- and geometrically-corrected, 
and dither-combined (via AstroDrizzle), the flc images 
have not been resampled to correct for geometric distortion 
and thus allow for more accurate PSF fitting. 

In brief, we first measure centroid positions for bright, 
relatively isolated stars in the drc images, since these im¬ 
ages have an astrometric solution that al lows us to con¬ 
struct an initial reference frame. As in lAnderson et al.l 
J2008al f. our reference frame has a 50 mas pixel 1 scale 
and a north-aligned y axis. Next, we perform PSF pho¬ 
tometry on the indi vidual flc exposures usin g the program 
img2xym_WFC. 09x10 l| Anderson fc Kingll200fil ). which uses a 
library of models of the spatially-variable effective ACS PSF. 
Once me asured, ste l lar po sitions are corrected for distortion 
using the lAndersonl (|2006i f corrections. Stars in common be¬ 
tween the distortion-corrected frame of each flc exposure 
and the reference frame are identified, allowing us to deter¬ 
mine six-parameter linear transformations from the former 
to the latter. We iterate on this process once, replacing the 
centroid positions that initially defined the reference frame 
with the average transformed reference-frame positions from 
the flc PSF photometry, then recomputing the linear trans¬ 
formation from each flc image to the reference frame. 

We relate both epochs to a single reference frame, so 
the measured positions in each epoch can be directly com¬ 
pared. Relative proper motions for point sources (including 
the YSOs adjacent to HH 900) are simply the difference 
in average reference frame position between 2005 and 2014. 
The positional errors in each dimension for each epoch are 
calculated by dividing the rms scatter of the individual mea¬ 
surements from each exposure, divided by the square root 


of the number of exposures (N=2-4) used to compute the 
average position. 

Our HST observations are such that the western bow 
shock of HH 900 was observed in a different orbit than the 
globule, inner jet, and eastern bow shock. Each orbit con¬ 
sists of three overlapping pointings, with two dithered ex¬ 
posures per pointing; however, the overlap in area between 
orbits is small. This requires that we construct a separate 
reference frame for each orbit. These reference frames are 
defined by the positions of the stars and are not tied to an 
absolute proper-motion zero point. Therefore, all motions 
are measured relative to the average motion of the stars in 
the image. Because the HH 900 western bow shock was ob¬ 
served in a separate frame, there may be a systematic offset 
between the proper motion of the western bow shock and 
the rest of HH 900. We expect this offset to be < 1 km s _ , 
but certainly no larger than the stellar velocity dispersion 
of the Tr 16 cluster, or < 5 km s -1 . 

To measure the proper motions of extended features 
in HH 900, we create stacked images by resampling the 
fl c ex posures into our reference frame (see lAnderson et ali 
l2008al . for details of the stacking algorithm used). We then 
select bright jet features that do not change significantly in 
morphology between the two epochs and measure their mo- 
tio n using the modified cross-correlation techn i que described 


tion using tnc modined cross-correlation technique described 
by Currie et alJ J 1 99 ( 1 ) ; lHartigan et al.l (120011 1 ; iMorse et al.l 
J200lt) . After subtracting a median-filtered image to account 
for background emission, we select a small box around each 
jet feature. This subimage is shifted relative to the refer¬ 
ence image, and for each shift, we compute the total of the 
square of the difference over the box region between the two 
images. The minimum of this array corresponds to the pixel 
offset between the two images. For the typical signal-to-noise 
in our images, the uncertainty of this procedure is ss 0.03 
pixel, corresponding to « 10 km s -1 . 

We also measure the proper motions of point sources in 
the image by fitting a spatially variable PSF to the stars. 
This allows the motion of stars to be measured to greater 
precision than nebulous jet features, typically a few km s -1 
(see Table 2). 

2.2 Spectroscopy 

FIRE: We obtained near-IR spectra of HH 900 with the 
Folded-Port InfraRed Echellette (FIRE) near-IR spectro- 
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Figure 1. New narrowband [Fe n] images of HH 900 obtained 
with HST /WFC3- IR. Panel (a) is the HST /ACS Ha image from 
Smith et al.l d2010l) showing the full extent of the optical out¬ 
flow. Panel (b) is a [Fe n] 1.26 pm + [Fe ii] 1 .64 pm image 
of HH 900 with the two YSOs identified by IShinn et alj (120131) 
from the PCYC labelled. Panels (c) and (d) show the continuum- 
subtracted [Fe ii] 1.26 pm and [Fe ii] 1.64 pm images, respectively. 


graph (|Simcoe et al.l 120081 , l2010l , l2013h on the Magellan 
Baade 6.5-m telescope on 17 January 2014 (see Table [l]). 
FIRE’s 0.8 — 2.5 pm wavelength coverage includes multiple 
emission lines from the jet. We positioned the 7" slit along 
bright [Fe ii] emission in the eastern and western limbs of 
the outflow and used a 1" slit width to accommodate > 11 
seeing. Figure [6] shows the slit positions. The YSO and pu¬ 
tative microjet that lie along the western limb of the larger 
bipolar outflow fall within the slit used to observe the west¬ 
ern limb of the [Fe n] jet. To account for sky emission, we 
employed a nodding strategy, pointing on- and off-source in 
an ABBA sequence with an average sky offset located ~ 2" 
below the jet. Wavelength calibration was done using the 
internal ThAr lamp and data reduction was performed with 



offsect [arcsec] 

Figure 2. Ha HST/ ACS image from [Smith et aD ([2010 ) (top) 
and narrowband H 2 — K-band image from GSAOI (bottom), 
both with [Fe 11 ] contours (1.26 pm + 1.64 pm) overplotted. [Fe 11 ] 
emission begins at the same place that H 2 emission fades below 
our sensitivity and the Ha emission narrows to a similar width 
as the [Fe 11 ] emission. 


the Firehose IDL pipeline. After correcting for the motion 
of the earth, we report Doppler velocities in the heliocentric 
velocity frame. 

EMMI: We obtained visual-wavelength spectra of 
HH 900 with the ESO Multi-Mode Instrument (EMMI) on 
the New Technology Telescope (NTT) on 09 March 2003. We 
observed HH 900 in echelle mode, using the cross-dispersing 
grism to obtain Ha and [S 11 ] simultaneously. The echelle slit 
length is 25” and the jet was observed with a 1" slit width. 
We aligned the slit through the middle of the globule and 
the star in the western limb of the flow. We estimate sky 
emission from separate frames, with the slit offset ~ 4” per¬ 
pendicular to the jet. We compute the wavelength solution 
from the sky frames, using the nebular [N 11 ] lines for the 
Ha dispersion solution and the [S 11 ] lines themselves. Thus 
our derived radial velocities are relative to the Doppler shift 
of ambient ionized gas in the Carina H 11 region. 


3 RESULTS 

3.1 Morphology in the IR images 

Figure [Q shows the new [Fe 11 ] 1.26 pm and 1.64 pm images 
of HH 900 obtained with F5T/WFC3-IR (Figures □Jo , c, and 
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Figure 3. Top: A [Fe n] 1.26 pm/ 1.64 pm ratio image (darker 
colours correspond to a higher value of IV) with lines indicating 
the boundaries used to make the [Fe n] ratio tracings along the 
length of the jet below plotted below. Middle: Tracing of the inner 
jet along the [Fe n] jet axis with positions plotted as a function 
of the offset from the approximate location of the center of the 
globule. Bottom: Same tracing as above, but with the two sides 
of the jet overplotted on top of each other. 


d) and the ACS Ha image (Figure |TJl). Bright, collimated 
[Fe n] emission that extends to the east and west of the 
dark tadpole globule stands out in the continuum-subtracted 
images (Figures [lj: and d). Near-IR [Fe n] emission traces a 
symmetric, collimated bipolar jet, unlike Ha, which emerges 
into the H n region with almost the same geometric width 
as the globule (see Figure [2j. Figure [2] shows a more detailed 
comparison of the Ha and [Fe n]. Intensity tracings through 
the eastern and western limbs of the [Fe n] jet show that 
they are remarkably similar (see Figure [3] and [5]). 


plot direction 

y-offset 
n r 

☆ 

_J 


-0.1" 



[Fe II] 1.26 pm /1.64 pm 



[Fe II] 1.26 pm /1.64 pm 


Figure 4. Top: [Fe n] 1.26 pm/ 1.64 pm ratio image with lines 
showing the boxes used to extract the tracings through the jet. A 
star in the middle denotes the estimated location of the invisible 
driving source. Bottom left: Plots of the average [Fe n] 1.26 pm 
/ 1.64 pm ratio through the width of the eastern limb of the jet 
(perpendicular to the direction of jet propagation). We sum the 
emission along the length of the jet, as shown in the image at top. 
Positive offsets on the y-axis are located above the jet, negative 
offsets fall below. Bottom right: The average [Fe n] 1.26 pm / 
1.64 pm ratio through the western limb of HH 900. 


Both sides of the jet show a > 1.5" gap between the 
edge of the globule and bright [Fe n] emission from the jet. 
Unlik e the images of jets in our previous study using archival 
data (jReiter Sz Smithll2013l F our observations of HH 900 in¬ 
clude simultaneous line-free continuum images that allow us 
to subtract the continuum and therefore separate photodis¬ 
sociation region (PDR) emission along the globule edge from 
fainter [Fe n] emission emitted by the inner jet. Subtract¬ 
ing the line-free continuum from HH 900 demonstrates that 
the offset of the [Fe n] jet from the edge of the globule is 
real (see Figure [2j). No tenuous [Fe n] emission from the in¬ 
ner jet extends back to the globule edge in the continuum- 
subtracted image, demonstrating that the offset is not just 
confusion with continuum emission that grows increasingly 
brig ht close to the irrad iated edge of the dusty globule. 

ISmith et al.1 (2010) identify two bow shocks associated 
with HH 900 that are ~ 10" beyond the end of the inner 
jet in Ha images. Both bow shocks are also bright in [Fe il] 
emission (see Figure [l|. 

Three dark, filamentary streaks in the western limb of 
HH 900 are seen in a bsorp tion in the Ha image (see Fig¬ 
ure [TJl) . ISmith et al.l ll2010li interpret these dark filaments 
as extinction due to dust entrained by the jet along the 
walls of the outflow ca vity, similar to the features that 
IWalawender et all (120041 ') observed in HH 280 (another jet 
that appears to burst out of a globule). Bright emission ex¬ 
tending off the western edge of the globule in the [Fe il] 
image in Figure [ljr disappears in the continuum-subtracted 
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image, indicating that it is continuum emission. This may 
be starlight scattered by the dusty filaments. 

Narrowband H 2 images from GSAOI show H 2 emission 
from the surface of the globule and outside it that, like the 
Ha emission, appears broad with nearly the same projected 
width as the globule itself (see Figure 0. The intensity of 
the H 2 emission appears to be slightly (~ 1.5x) greater 
along the edges of the outflow than in the center (i.e. limb 
brightened). H 2 emission from the outflow fades below our 
detection limit ~ T/5 away from the globule edge; this is the 
same place where the width of the Ha emission appears to 
decrease and the collimated [Fe 11 ] jet begins (see Figure [5]). 

3.2 [Fe 11 ] ratio 

Both the [Fe 11 ] 1.26 pm and 1.64 pm lines originate from 
the same upper level, so the observed flux ratio therefore 
provides a measure of the reddening. In unobscured envi¬ 
ronments, the 1.26 pm line will be brigh ter, wi th a derived 
flux ratio, TZ = A12567/A16435 = 1.49 llSmith fe Haiti aarl 
I 2 OO 6 I I. Smaller values of TZ indicate more reddening and ex¬ 
tinction. Along the length of the [Fe 11 ] jet, the ratio remains 
consistent, with a value of ~ 0 .95 (see Figured). Th is corre¬ 
sponds to an E(B—V) « 0.9 (iReiter fc Smithj[2013li . Using 
the ratio of total t o selective extinction measur ed toward 
Carina of R = 4.8 (|R. Smithlll987l : ISmithll2002T ). this cor¬ 
responds to Ay ~ 4.5 mag along the length of the jet. We 
also measure TZ through the width of the jet (perpendicular 
to the direction of propagation, see Figure HI). Reddening 
through the width of the jet increases further away from rj 
Car and Tr 16, corresponding to a decrease in TZ of ~ 0.2 
across the width of the jet. The measured decrease in TZ cor¬ 
responds to ~ 1.5 mag more extinction on the side of the jet 
further away from the ionizing source. 

3.3 Candidate driving sources 

[Reiter fe Smithl ([20131 1 showed that in HH 666 and HH 1066, 
the driving protostar lies along the projected jet axis and 
that [Fe 11 ] connects the jet to the Spitzer -identified driving 
sources. However, [Fe n] emission in HH 900 does not reach 
back to the edge of the tadpole-shaped globule, and no [Fe 11 ] 
emission from the jet is detected inside the globule. 

Two IR-bright YSOs h ave be en found in the vicinity of 
the HH 900 globule dShinn et alj|2013t) . The jet axis traced 
by [Fe n] emission lies ~ 1.5" north of the YSO that lie s at 
the bottom of the tadpole globule (PCYC 842: lPovich et all 
l201ll. see Figure [!}, clearly demonstrating that it cannot be 
the driving source of the bipolar flow. The second Spitzer- 
identified YSO (PCYC 838) lies along the western limb of 
the flow, in the middle of the putative microjet. While the 
morphology raises the possibility that this protostar drives 
a separate jet that happens to be aligned with the larger Ha 
outflow lobe, images alone do not indicate whether there is 
a separate microjet or if PCYC 838 is its driving source. 
As we discuss in Section 13.51 however, radial velocities and 
proper motions lead us to reject this hypothesis. 

Evidence for an additional protostellar source along the 
jet axis embedded in the tadpole globule remains elusive due 
to the small size of the globule (~ 2") and the coarser reso¬ 
lution of available mid- and far-IR observations. We do not 




offset from globule center [arcsec] 

Figure 5. Top: Median intensity traced along the axis of the 
[Fe 11 ] jet in HH 900 (using the same tracing area shown in Fig¬ 
ure Bottom: Median intensity in the eastern and western limbs 
of the jet over plotted on each other. 
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detect any evidence of an IR-bright point source inside the 
globule in any of the WFC3-IR images. The upper limit in 
all four bands is ^ 21 mag. High frequency, high spatial res¬ 
olution observations with ALMA are required to determine 
whether there is an additional protostar embedded in the 
globule. 


3.4 Velocity structure 

Among the peculiarities of HH 900 revealed with Ha images 
from HST is the morphology of the dark globule. A small, 
wiggly, dark tail extends from the southeastern edge of the 
globule, creating a tadpole-like morphology (see Figure |T|). 
The tail of the dark tadpole globule appears to lie in front of 
the jet, leading iSmith et alj (|20ld l to suggest that the east¬ 
ern limb of the jet is redshifted, and the western limb of the 
jet is blueshifted. Both optical and IR spectra now confirm 
this conjecture (see Figures [6] and [7]). Ha emission from the 
western limb of the jet is blueshifted and shows a Hubble- 
like velocity structure, with Doppler velocities extending up 
to roughly —40 km s^ 1 . A similar velocity structure is barely 
discernible in the [S il] spectrum (the average of the A6717 
and A6731 emission) due to the lower signal-to-noise. 

In contrast to the Hubble-like flow in the optical emis¬ 
sion, near-IR. [Fe il] lines from both sides of the jet show a 
relatively constant velocity with redshifted velocities mea¬ 
sured to be ~ 30 — 60 km s -1 in the eastern limb of the jet 
and blueshifted velocities of ~ 0 — 20 km s -1 in the western 
limb (suggesting v sys « 19 km s -1 , see Figure[6]). The veloc¬ 
ity difference between the eastern and western limbs of the 
jet is symmetric about the globule, indicating that the [Fe il] 
emission is dominated by a single bipolar jet launched from 
within the globule. The H 2 emission from the eastern and 
western limbs of the jet is also redshifted and blueshifted, re¬ 
spectively. However, H 2 emission appears to be ~ 20 km s -1 
bluer than [Fe 11 ] velocities in both limbs. Position-velocity 
diagrams illustrate that the locations of H 2 and [Fe 11 ] emis¬ 
sion from the jet appear to be mutually exclusive (see Fig- 
ure[S]). In fact, Ha, [Fe 11 ], and H 2 emission all trace different 
velocity components (see Figures [6] and [7]). Blueshifted ve¬ 
locities in the Ha position-velocity diagram increase from 
~ 0 km s -1 to the higher velocity traced by the [Fe 11 ] jet 
(~ —40 km s -1 ) while H 2 emission appears to be ~ 20 km 
s -1 bluer than v ays inferred from the [Fe 11 ] emission from 
the innermost part of the jet. 


3.5 Ha proper motions and 3-D velocities of jet 
features 

With a ~ 9 yr time baseline between the two HST /ACS 
epochs, we are sensitive to transverse velocities of ~ 25 km 
s' 1 in nebulous jet features seen in Ha, limited by the signal- 
to-noise of the jet knot compared to the bright background 
of the H 11 region. We report proper motions and radial ve¬ 
locities in Tabic [2] Despite the relatively long time baseline, 
the morphology of the outflow has remained remarkably con¬ 
stant. Only the bow shocks and knots capping either side of 
the broad bipolar outflow show noticeable motion between 
the two epochs. Both knots move away from the globule 
along the jet axis with transverse velocities ~ 100 km s -1 
(see Figure ©. The Ha-bright filament along the western 
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limb of the jet does not appear to move at all (< 25 km 



The eastern and western bow shocks move in opposite 
directions, with proper motion velocities near the tip of the 
shock of ~ 60 km s -1 . For jets that lie near the plane of the 
sky (as appears to be the case for HH 900, see Table [2]), the 
difference between the inclination-corrected proper motions 
of adjacent knots can be used to estimate the shock veloc¬ 
ity. Both bow shocks have velocities ~ 40 km s -1 slower 
than the inner jet, suggesting shock velocities ~ 40 km s” 1 . 
Fast, dissociative J-shocks ( v s hock > 30 — 40 km s -1 ) a re ex¬ 
pected to be brigh t in near-IR [Fe il] lines llNisini et al .1120021 : 
IPodio et al .1120061 ) and indeed, both shocks have bright [Fe 11 ] 
emission (see Section E and Figure [lj. 

Combining the transverse velocity of jet knots with the 
radial velocity from spectra, we can calculate the tilt angle, 
a = tan ~ 1 (v r /vT) and 3-D space velocity. With the low 
radial velocities in the Ha spectrum, we find that the jet 
lies close to the plane of the sky (although we note the poor 
signal-to-noise of the spectrum). We derive a tilt angle of the 
jet away from the plane of the sky of < 10° (see Table ®. 


3.6 Protostar kinematics 

We can also constrain the motion of the two protostars near 
HH 900. The YSO that lies along the western limb of the 
outflow (PCYC 838, see Figure [1} falls within the slit we 
use to observe the western limb of the jet. We extract the 
spectrum of this source separately. After subtracting off ex¬ 
tended emission, we find that the central velocity of the 
hydrogen recombination lines in the spectrum of the YSO 
constrain the heliocentric Doppler velocity of the protostar 
to be < |5| km s _1 (compared to the heliocent ric sy stemic 
velocity of —8.1 ± 1 km s ^ 1 found for 77 Car bv lSmithll2004| j. 

We also measure the proper motions of both protostars. 
PCYC 838 moves in roughly the same direction as the jet 
with a transverse velocity of 4.3 ± 1.3 km s _1 (see Table 2). 
Together, the radial and transverse velocities constrain the 
motion of PCYC 838 to < 7 km s -1 . We do not have spec¬ 
tra of the protostar near the bottom of the dark globule, 
PCYC 842 (see Figure [I}, but we can measure its proper 
motion in the plane of the sky. The projected motion of 
PCYC 842 is toward the globule with a transverse velocity 
of 15.4 ± 3.0 km s -1 (see Figure [HJ. 


4 DISCUSSION 

4.1 Spatial offset of [Fe 11 ] emission 

Bright [Fe II] emission on either side of the dusty tadpole 
globule in HH 900 traces a symmetric, narrow jet. Like 
HH 901 and HH 902, two externally irradiated jets that are 
also in the Carina nebula, the brightest [Fe 11 ] emission turns 
on more than 1" away from the edge of the globule, indi¬ 
cating a gap between the globule edge and the beginning of 
bright [Fe 11 ] emission from the jet. The offset of the [Fe 11 ] 
emission may be interpreted as evidence that the mass-loss 
rate in the jet has recently decreased, since a lower den¬ 
sity jet may not maintain a column of material sufficient t o 
shield Fe + from further ionization llReiter fe Smith| l2013l b 
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Figure 6. From left to right: FIRE slit positions drawn on a [Fe ii] 1.64 image, alongside position-velocity diagrams of HH 900 in the 
[Fe n] 1.26 pm, [Fe n] 1.64 pm, and H 2 S(0) 1-0 2.12 pm lines, respectively in the heliocentric velocity frame. Plots on the far right show 
the median line velocity over steps of ~ 0 / /6. Dashed lines indicate the systemic velocity inferred from velocity difference in the [Fe 11 ] 
emission from the eastern and western limbs of the jet. Continuum emission in spectra of the western limb of the jet (top row) is from 
the protostar in the middle of the putative Hq microjet. [Fe 11 ] emission from the jet appears to be remarkably symmetrical on either 
side of the continuum. PV diagrams illustrate that the eastern and western limbs of the jet are redshifted and blueshifted (respectively) 
and the [Fe 11 ] emission grows bright where the H 2 emission fades away. 
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Figure 8. New HST /ACS Ha image of HH 900 with boxes used 
to measure proper motions overplotted. 


However, HH 900, HH 901, and HH 902 all have gaps of sim¬ 
ilar size, and a nearly synchronized end of an outflow phase 


from three physically unrelated jets would seem highly sus¬ 
picious. The greatest similarity between the three sources is 
their proximity (< 3 pc in projection) to a young massive 
cluster with log(Qn)> 50 photons s _1 (Tr 14 for HH 901 
and HH 902, Tr 16 for HH 900). This argues that the gap 
is produced by something in the environment, rather than 
due to a recent change in the mass-loss rate. 

The thre e dark , narrow str eaks in the western limb of 
HH 900 that ISmith et al] (|2010l j interpret as dust filaments 
in the walls of the outflow cavity suggest another explana¬ 
tion for the offset [Fe n] emission. These dusty streaks hint 
at a potentially large column of cold, dusty material that 
was blasted out of the globule by the jet. Close to the cloud 
edge, recently ejected dust and molecules may prevent UV 
emission from illuminating the inner jet. 

We can measure reddening in the HH 900 jet using the 
ratio TZ = A12567/A16435. 77 increases through the width 
of the jet (perpendicular to the direction of propagation), 
supporting the idea that external irradiation from massive 
stars in Tr 16 dominate the ionization and dust survival in 
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Figure 7. Ha and [S II] position-velocity diagrams of HH 900 shown on either side of an Ha image with the EMMI slit position show 
at approximately the same spatial resolution. Plots show the median velocity taken over ~ 2" increments. 


Table 2. HH 900 Proper Motions 


Object 

<5x 

mas 

mas 

V T a 

[km s -1 ] 

Vi? b 
[km s -1 ] 

velocity 0 
[km s” 1 ] 

a 

[degrees] 

age d 

yr 

HH 900 A 

43 (i) 

-24 (0.4) 

60 (1) 


60 (1) 


3523 (115) 

HH 900 B 

31 (1) 

-2 (2) 

37 (2) 


37 (2) 


5364 (307) 

HH 900 C 

55 (0.1) 

-11 (1) 

68 (1) 


68 (2) 


2803 (87) 

HH 900 D 

74 (1) 

-31 (4) 

97 (3) 


97 (3) 


859 (34) 

HH 900 E 

76 (2) 

-25 (1) 

97 (3) 


97 (3) 


682 (25) 

HH 900 F 

-60 (1) 

42 (2) 

89 (2) 

-16 (4) 

90 (5) 

-10 (0.3) 

691 (24) 

HH 900 G 

-77 (1) 

39 (1) 

104 (2) 

-7 (10) 

105 (10) 

-4 (0.1) 

746 (24) 

HH 900 H 

-40 (3) 

14 (1) 

51 (3) 


51 (3) 


5379 (365) 

HH 900 I 

-42 (1) 

28 (1) 

61 (2) 


61 (2) 


4664 (174) 

PCYC 838 

-3 (1) 

2 (1) 

4 (1) 


4 (1) 


11201 (3126) 

PCYC 842 

5 (3) 

-11.8 (0.4) 

15 (2) 


15 (2) 


1393 (140) 


Proper motions measured for the YSOs and the knots marked in FigureJS] 

Uncertainties are listed in parentheses alongside the best-fit value. 

“ The transverse velocity, assuming a distance of 2.3 kpc. 
b The radial velocity measured from spectra. 

c The total velocity, assuming the average inclination when a radial velocity is not measured directly. 
r! Time for the object to reach its current location at the measured velocity, assuming ballistic motion. 


the jet, not the nearby protostars or shocks. The increase in 
1Z also suggests that some entrained dust survives even in 
the outer regions of the jet. 

Dust entrained in the outflow may obscure the inner jet, 
explaining the offset between [Fe n] emission and the edge 
of the globule. We propose that [Fe il] emission in HH 900 
traces the collimated jet while broad Ha and H 2 emission 
near the globule trace the wider body of an externally illumi¬ 
nated outflow cavity (see Figure [9]). Narrow [Fe 11 ] emission 
tracing the body of the jet has been s een in other HH jets i n 
Carina (e.g. HH 666 and HH 1066, see lReiter fc Smithl2013l) . 
Broad H 2 emission in HH 900 extends beyond the edge of 
the globule along the outflow direction, but only to the point 
where [Fe 11 ] emission begins, indicating that molecules per¬ 
sist in this intermediate zone. Extinction in the Ha image 
also hints at the presence of entrained dust that survives in 


the H 11 region close to the globule. Continuum-subtracted 
[Fe 11 ] and H 2 images are required to determine whether en¬ 
trained dust and molecules also obscure the inner regions of 
HH 901 and HH 902. 

HH 900 has two bow shocks separated from continuous 
emission in the inner jet, demonstrating that there has been 
a previous breakout of the jet into the H 11 region. Despite 
this, new near-IR images show evidence for a significant col¬ 
umn of dust and molecules outside the globule entrained 
by the jet. This may be expected from the youngest sources 
that are still surrounded by a substantial m olecular envelope 
with only a small cavity cleared by the jet (lArce fe Sargentl 
1200611 . Significant entrainment of molecular and dusty mate¬ 
rial by jets driven by more evolved protostars, or after many 
jet bursts, may be a direct consequence of the feedback- 
dominated environment (see Section 14.611 . Compression of 
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the molecular globule by massive star feedback may refill 
the cavity along the jet axis, increasing the amount of ma¬ 
terial entrained by a subsequent pulse of the jet. 

Both HH 901 and HH 902 also show evidence for previ¬ 
ous jet activity. If entrained material obscures the inner jet 
in all three cases, then this argues for a large column of mate¬ 
rial t o be drag ge d out of the globule repeatedly by jet bursts. 
iReiter fe Smithl d2013h estimate uh > 10 s cm -3 for the pil¬ 
lar housing HH 901, and we expect similarly high densities 
for the HH 902 and HH 900 clouds (also see Section 14.411 . 
All three jets are embedded in the hot stellar wind bubble 
created by the many O-type stars in the nearby Tr 14/Tr 16 
star cluster, so they are subject to similar amounts of feed- 
ba ck. Jets emerging from more diffuse globul es (e.g. HH 666, 
see lSmith et al.lT2004al : IReiter fe Smithfl2013l . Reiter et al. in 
preparation) do not show this same offset from the globule 
edge. 

4.2 The putative Ha microjet 

Based on the morphology in the Ha image, ISmith et al.l 
d20ld! identify a possible microjet along the western limb 
of the broad HH 900 outflow. A bright point source with 
colours consistent with being a YSO lies in the middle of 
this bright Ha filament and was identified as the possible 
driving source for the putative microjet. After examining 
the new imaging and spectroscopy presented here, we find it 
unlikely that this Ha-bright feature is a separate microjet. 
The key arguments against the microjet hypothesis are (1) 
the morphology in the [Fe n] images, (2) radial velocities in 
spectra, and (3) proper motions as detailed below. 

(1) The [Fe n] emission indicates that the main body of 
the HH 900 jet is symmetric about the center of the globule 
and shows no significant deviation at the position of the pu¬ 
tative microjet. [Fe n] intensity tracings through both sides 
of the jet are similar, with no significant difference in the 
brightness of the western limb as compared to the east (see 
Figure [3] and Section [3.111 . It is possible that the microjet, 
driven by an unobscured star, is completely ionized, and 
therefore is only visible in Ha. However, if this were the case, 
the jet would have to be projected in front of the blueshifted 
limb of HH 900. 

(2) Linewidths of the [Fe il] and H 2 emission from the 
western limb are similar to the eastern side of the jet, and 
appear to trace a single jet (see Figure [S] and Section [3.411 . 
Moreover, velocity changes on either side of the jet are mono¬ 
tonic. There is no evidence for two jet velocity components 
in the Ha spectrum (although the signal-to-noise and veloc¬ 
ities are low, see Figure 0. 

(3) The high degree of symmetry in the Ha filament 
on either side of the protostar combined with the lack of 
any radial velocity changes suggests that if this feature is 
indeed a jet, it must lie almost exactly in the plane of the 
sky. However, our proper motion measurements constrain 
any oppositely directed, outward motions in the plane of the 
sky to < 25 km s _ . Together, proper motions and spectra 
require that any outflowing gas must move slower than ~ 30 
km s“ , unusually slow for a jet driven by an unobscured 
proto star (see, e.g. iReipurth fe Ballvll200ll ; IReiter fe Smithl 
l2014lh 

Without evidence for any bipolar motion in the Ha fil¬ 
ament away from PCYC 838, it seems unlikely that this 


putative “microjet” is a separate jet at all. Hypothetically, 
the YSO that appears to lie at its center may have been 
ejected from the globule in a dynamical encounter, with the 
motion of the ejected star through the larger bipolar out¬ 
flow creating a bright Ha tidal tail that looks like an Ha 
microjet. However, hydrogen emission lines in the spectrum 
of the YSO (seen intersecting the jet in Figure 0 constrain 
the Doppler velocity of the protostar to be v sya < —5 km 
s -1 whereas the Ha emission at the same position is ~ — 12 
km s _1 . Moreover, we can constrain the proper motion of 
the YSO with the two epochs of HST / ACS data. We find 
that the star moves almost imperceptibly in the Ha im¬ 
ages, with a measured velocity in the plane of the sky of 

4.3 ± 1.3 km s _1 . The YSO is not moving faster than the 
random motion of stars in the Carina nebula, making the 
ejection scenario unlikely. At these low velocities, the dy¬ 
namical time required for the star to reach its current po¬ 
sition is ~ 11,000 yr, longer than the estimated dynamical 
age of the two outer bow shocks in HH 900 (~ 2200 yr, see 
Section 13.Ill , and much longer than the age of the inner jet 
(~ 1000 yr, see Section n~n . Given the discrepancy in the 
timescales for the star and the jet, any interaction between 
the star and HH 900 is likely due to the jet recently mov¬ 
ing past the position of the YSO. The simplest explanation 
is a chance alignment of the YSO with the apparent center 
of the Ha-bright filament. Indeed. [Wolk et ~ahl (|201ll f find a 
relatively high spatially averaged source density in the Tr 16 
subcluster near HH 900 (45 sre pc -2 in sub-cluster 12). 

The Ha-bright filament may result from local irradia¬ 
tion due to the small separation between the YSO and the 
HH 900 outflow. Local irradiation from the driving protostar 
has been invoked to explain bright optical and IR emission 
from the inn er regions of protoste llar jets in more quiescent 
regions (e.g. iReipurth et aHl2000l '). The IR-bright YSO that 
appears to lie in the middle of this filament must lie ex¬ 
tremely close to the broad, bipolar outflow (< 4.5 x 10~ 4 
pc or < 92 au in the plane of the sky). Uncertainties in 
the three-dimensional geometry of the system may hide a 
slightly greater separation, although the presence of espe¬ 
cially bright Ha emission on either side of the YSO suggests 
that the protostar may affect the larger, apparently separate 
outflow. Thus, while our kinematic data rule out the possi¬ 
bility that this bright Ha feature is a separate microjet, its 
origin remains elusive. 

4.3 Mass-loss rate in the neutral jet 

[Smith et al.l (j2010l ~) estimated mass-loss rates, Mjet, of the 
jets in Carina using the Ha emission measure in ACS im¬ 
ages. For HH 900, the wide inner jet in the east and the Ha- 
bright filament (the now refuted microjet) in the west repre¬ 
sent the two highest Mj e t in the sample. However, mass-loss 
rates calculated from the Ha emission measure reflect the 
total Mj e t only for fully ioniz ed jets , and otherwise provide 
a lower-limit. IReiter fc Smithl (120131 1 argue that the HH jets 
in Carina are not fully ionized, and that [Fe 11 ] emission 
traces neutral material in jets that are sufficiently dense 
to self-shield (i.e. to prevent Fe + —> Fe ++ ). Relaxing the 
assumption that the jets are fully ionized and calculating 
the density required to self-shield generally increases the es- 
timated mass-loss rat e by at least an order of magnitude 
llReiter fe Smithll2013ll . For HH 900, the situation is unfor- 
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tunately even more complex. Ha and [Fe n] have different 
spatial distributions and appear to trace different velocity 
components, suggesting that the mass-loss rate estimated 
from Ha only samples the mass contained in a thin ionized 
skin on the outside of the outflow sheath entrained by the 
[Fe n] jet. This is remarkable since the mass-loss rate derived 
from the Ha emissi on measure was al ready the highest of all 
the jets in Carina rtSmith et al.ll2ofol l. 

We can estimate the mass-loss rate of the HH 900 
[Fe n] jet by requiring that a neutral, cylindrical jet survives 
out to a distance Li before it is comp letely photoablated 
llBallv et al.ll2006l : iReiter fe Smithll2013l f . Here, we take L\ 
to be the length of the continuous [Fe n] jet (not including 
the bow shocks). For a cylindrical jet, 


Mjet 


Liffim H cn 

2D 


OLB 

nrj et QHsin{/3) 


- 1/2 


(1) 


where / ~ 1 is the filling factor, fi ss 1.35 is the mean molec¬ 
ular weight, niH is the mass of hydrogen, cii ~ 11 km s -1 is 
the sound speed in ionized gas, and ag « 2.6x 10 -13 cm 3 s~ x 
is the Case B recombination coefficient. CPD-59 2641 is an 
05 V star located ~ 1' from HH 900, or a projected distance 
D ~ 0.7 pc (which may underestimate the true separation 
by ~ y/2). Assuming that CPD-59 2641 dominates the ion¬ 
ization, the i oniz ing photon luminosity is log(Qn) = 49.22 
s _1 (ISmithlfiood l. We measure the jet radius, Vj e t « 0.007 
pc, from the [Fe n] images and assume that the angle /3 be¬ 
tween the jet axis and the direction of the ionizing radiation, 
is 0 ~ 90°. From this, we find Mj e t > 7 x 1CU 6 Mq yr -1 
and Mj e t >5x 1CU 6 Mq yr -1 , respectively, for the eastern 
and western limbs of the jet (compared to 5.68 x 10“' Mq 
yr -1 and 6.20 x 10“' Mq yr -1 , respectivel y , esti mated from 
the Ha emission measure. seelSmith et ahll2010l l. 

As IReiter &; Smithl d2013f ) found for other HH jets in 
Carina, using [Fe n] as a diagnostic and without assum¬ 
ing the jet is fully ionized, M j e t is an order of magnitude 
larger than lSmith et all (120101 1 estimate from the Ha emis¬ 
sion measure. Assuming that the mass accretion rate, M acc , 
onto the protostar is ~ 10 — 100 times larger t han the mass - 
loss rate in the jet (e.g. IllaiTiaan et al.lfl995l : ICalved[l99Sll . 
we find M acc ~ 10 -4 Mq yr -1 . This points to a high lumi¬ 
nosity from the driving protostar, either due to its relatively 
high mass or a recent accretion burst. 


4.4 Protostars and globule survival 

The jet axis defined by the collimated [Fe n] emission rules 
out the hypothesis that the driving source of the jet is the 
pr otostar at th e bottom of the HH 900 globule (identified 
by ! Shinn et iDlmi). Although the jet axis dehned by the 
[Fe n] emission in the WFC3-IR images from HST indicates 
a clear mis-match of the star and the jet (see Figure [2]), is it 
possible that the star began in the globule and was ejected, 
perhaps in a small-N stellar inte raction fe.g. lReipurth et al.1 
l20ld : fReipurth fe Mikkolall2012l l? If this is the case, the star 
must have been in the globule at the time that the innermost 
jet emission we detect in [Fe n] was ejected. Assuming a 
jet velocity of 100 km s _1 for the [Fe ii] emission located 
~ l"5 from the globule edge (thus ~ 2'/5 from the presumed 
position of the driving source), this means the star would 
need to have been ejected within the last ~ 275 yr. To travel 
to its current location at the bottom of the globule in that 


time, the star would have to be moving at > 60 km s _1 . 
However, PCYC 842 has an observed proper motion of 15.4 
km s" 1 toward the globule (see Section 13.61 and Tabic [2]), 
clearly inconsistent with the ejection hypothesis. 

Without compelling evidence for the recent ejection of 
the jet-driving source, we must conclude that an additional 
protostar remains hidden inside the dark globule. Symmet¬ 
ric [Fe ii] emission from the eastern and western sides of 
the jet, including a similar offset from the globule edge (see 
Figure 0, argues that the driving source is located inside 
the globule. Dark lanes in the blueshifted western limb of 
the jet apparent in the Ha image connect the jet directly 
with the globule and require that it is driven by a protostar 
embedded within it. However, we do not detect IR emission 
associated with a protostar in the globule. 

The high mass-loss rate in HH 900 ma kes t he non ¬ 
detect ion of a protostar particula rly puzzling. ISmith et al.1 
(l2010h and IReiter fe Smithl (l2013l l have argued that the HH 
jets in Carina are driven by intermediate-mass (~ 2 — 8 
Mq) protostars based on their high mass-loss rates, and 
the luminosity of candidate driving sou rces f ound fo r some 
of the jets support th i s interpretation (Smith et_al.l l2004al : 
lOhlendorf et alJ 120121 : IReiter fe Smithl 20131 . 2014ll . Even 
if the HH 900 driving source were a low-mass protostar, 
it would need to be undergoing an FU-Orionis-type out¬ 
burst to create a j et of this strength (ICroswell et al.l[l987l : 
ICalvet et al.lll993f) . and would therefore likely have a high 
accretion luminosity. However, we note that the dynamical 
age of the inner jet is ~ 1000 years (for a jet velocity of 
100 km s _1 ), a factor of ~ 10 longer th an t h e ty pical de¬ 
cay t ime of an FU Orionis outburst dHartmann fc Kenvonl 
11996h . Small velocity differences between the redshifted and 
blueshifted sides of the jet suggest that the outflow axis lies 
close to the plane of the sky, and we find a small inclina¬ 
tion angle for the jet (a < 10°). Thus the circumstellar disc 
of the driving source will be seen almost edge-on. A YSO 
viewed through the midplane of an optically thick disc will 
be heavily extincted by a large column of circumstellar ma¬ 
terial. If this YSO is embedded in a dense and sufficiently 
opaque globule, no scattered light will escape the cloud and 
the protostar will remain unseen. 

To estimate the amount of circumstellar material re¬ 
quired to completel y obscure the dr i ving s ource, we use the 
YSO models from iRobitaille et al.l (l2006ff . Since we have 
only upper limits on the protostar postulated to be in the 
globul e (see Section m, we instead do a parameter search 
of the IRobitaille et al.l (120061 ) models. We make the conser¬ 
vative estimate that the protostellar mass is 1.95 — 2.05 Mq 
and the opening angle of the envelope cavity angle is 5—10°. 
We estimate a small cavity opening angle based on the fact 
that Ha and H 2 emission are not wider than the globule. 
Within this search criteria, we sample disc accretion rates 
as high as ~ 10 -5 Mq yr -1 , on the lower end of the ac¬ 
cretion rates we expect if Mj e t ~ 0.01 — 0.1 x M aC c (see 
Section S3J. The disc accretion rate will be even higher 
(leading to a higher accretion luminosity) if the driving 
source is more massive (closer to 8 Mg), or undergoing an 
FU Orion is outburst (with disc accretion rates > 10 -4 Mq 
yr -1 , e.g. lHartmann fe Calvetlll995l ; ICalvetl 119981 1. For the 
12 models that satisfy this search criteria, the estimated 
line-of-sight extinction to the stellar surface ranges from 
A v fs 5 x 10 3 — 2 x 10 7 . Allowing for larger cavity opening 
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angles of 40 — 45° corresponding to more evolved driving 
sources, A v may be as low as « 3 x 10 2 . Taking the lower 
bound, A v = 10 2 , we can make a crude estim ate of the glob - 
ule mass. Using Nh ~ 1.8 x 10 21 x A v cm -2 (iGiiver &i O zell 
120091 1. the column density along the line of sight must be 
at least ~ 2 x 10 23 cm -2 . Assuming the globule is a uni¬ 
form sphere, we can estimate the volume density. The di¬ 
ameter of the globule is ~ 2", corresponding to a linear size 
of ~ 0.02 pc at the distance of Carina. For a globule ra¬ 
dius of 0.01 pc or 3 x 10 16 cm, this rough estimate yields 
an average volume density n ~ 6 x 10 6 cm -3 , among the 
higher volume densities inferred from dust emi ssion of low- 
mass star f orming cores (~ 10 4 — 10 6 cm -3 , lEnoch et al.l 
120071 . 120081 . a lthough higher den sity cores have also been 
observed, e.g. iTokuda et al.ll2014f ). Taking the globule to be 
a uniform sphere, this density suggests a globule mass of > 1 


M 0 . _ 

|Grenman fc Gahml (I2014T ) use the Ha images of 
I Smith et al.l 1 2010h toestimate the extinction in the HH 900 
globule and calculate a mass of ~ 14 Mj up iter, two or¬ 
ders of magnitude smaller than our estimate. This single¬ 
wavelength assessment of the extinction provides only a 
lower limit on the true column of material in the globule and 
is compromised by emission from the surface of the limb- 
brightened globule. Indeed, requiring that the globule has 
survived photoevapo ration f rom the many O-type stars in 
Tr 16 for > 3.5 Myr llSmithll2006l l argues for a much larger 
mass reservoir. The photoevaporation mass-loss rate of the 
globule is 


M phot « 47T r 2 n H ^m H cu 


( 2 ) 


where r is the radius of the globule, and tih is the density 
of neutral hydrogen. We can estimate the density at the 
ionization front from the requirement that ionizations are 
balanced by recombinations. This gives 

UIF “ \/ 8 nD?ra B (3) 

where, as in Section m we assume that the nearby O- 
type star CPD-59 2641 dominates the ionization. This yields 
niF ~ 4000 cm' 3 . We can connect this to the density in 
the globule using the continuity equation, tlifCii = uhCi 
where ci ~ 3 km s _1 is the sound speed in the molecular 
globule. For the estimated uh ~ 1.5 x 10 4 cm -3 , we find 
Mphot ~ 7 x I0~ 6 M@ yr~ 4 . This is within a fa ctor of two of 
the photoevaporation rate lSmith et al.l d2004bl f calculate for 
the nearby finger globule in Carina (~ 2 x 10“ 5 M 0 yr _1 ) as 
the higher flux of the ionizing source (log(Qrr) = 50.0 s _1 ) 
is roughly canceled by its larger separation (~ 6'). With this 
high photoevaporation rate, a 14 M j up globule will rapidly 
be ablated, and completely destroyed in ~ 2000 years. In 
contrast, the remaining lifetime of a > 1 M 0 globule will be 
> 10 5 years. 


4.5 Comparison to molecular outflows 

iReiter fc Smitbl (l20I3l j postulated that the irradiated jets in 
Carina are similar to the jets that drive molecular outflows, 
but are laid bare by the harsh UV radiation permeating Ca¬ 
rina. In this picture, molecules entrained in the outflow that 
are swept into the H 11 region will be quickly dissociated, 
and should only be observed within a few arcsec of the edge 



m 

wilds \ ^ ^ ^ 


Observational diagnostics 


bow shock from 
previous jet burst 


Figure 9. A time-series illustration of an outflow breaking out of 
a small, dense globule into the H II region. Panel I shows the jet 
turning on and beginning to entrain an outflow before emerging 
from the globule as shown in panel II. Panel III illustrates the 
current geometry of the irradiated jet-outflow structure. Different 
materials will dominate the observed structure of the outflow as it 
propagates further from the globule and spends more time subject 
to external irradiation (left side). Observational tracers of the 
different components of the irradiated outflow (right side) reveal 
the decreasing obscuration of the jet and increasing impact of 
external irradiation. Red lines trace [Fe ii] emission, blue lines 
show Ha, and gray lines illustrate the location of molecules. 


of the natal globule. Indeed, H 2 emission in HH 900 survives 
only ~ 1.5" away from the edge of the molecular globule, 
and must have been introduced into the H 11 region recently. 
Other authors have invoked a two-component outflow model 
to explain the coexistence of a high velocity jet and a slower, 
broader molecular outflow created by a wide-angle wind 
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(e.g. HH 111, N aear et al.l Il997l : iLee et ah 200d: H H 315, 
lArce &; Goodman 120021 : HH 46/47. lArce et all 20131 '). 

In HH 900, there is the further complication that the 
Ha position-velocity diagram traces a smooth increase in 
velocity with increasing distance from the dark globule. In 
the case of HH 111, higher velocity emission far from the 
driving source has been seen in maps of the CO emission 
and interpreted as evidence of jet-bow- shock entrainment 
(e.g. ILee et all 120001 : iLefloch et al.l l2007tl . As a bow shock 
propagates through its parent cloud, it will sweep up ma¬ 
terial, creating a sheath that envelopes the jet traveling in 
its wake. The initial impulse of the jet turning on (or up) 
creates the Hubble-like velocity structure. A corresponding 
shell of material swept up in HH 900 will be externally ir¬ 
radiated by Tr 16 upon exiting the globule, creating an ion¬ 
ized outer layer that will emit Ha. Within ~ 1.5" of the 
globule, this ionized sheath still exhibits H 2 emission. The 
overall blueshift of the H 2 velocities (~ 20 km s” 1 ) com¬ 
pared to [Fe 11 ] suggests that H 2 emission is dominated by 
the expansion of the outflow cavity. The abrupt end of H 2 
emission where [Fe 11 ] emission begins may signal a transi¬ 
tion in the column of obscuring, entrained material from a 
mix of atoms and molecules to primarily atomic gas. Beyond 
this point, Ha emission narrows until it converges with the 
[Fe 11 ] morphology at the head of the inner jet, as expected 
for prompt entrainment. Indeed, the Hubble-like velocities 
in the Ha position-velocity diagram suggest that HH 900 is 
the irradiated analog of an entrained molecular outflow. 


4.6 Environmental interaction 


Theoretical work on small globules in H 11 regions has 
explored their photoevaporation and the possibility that 
radiative compression may trigger star formation (e.g. 
Bertoldi 19891: Bertoldi fe McKee 1990l: Bertoldi fe Drainel 


19961 ; iGorti fe Hollenbachl 120021 ; iDale et al.l l2007allbl) . The 
small size of HH 900 and its survival despite its proxim¬ 
ity to rj Car hint that the protostar forming within it may 
have been triggered by radiation-driven implosion. There 
remains an ambiguity between the triggering of new star 
formation and the uncovering, and/or possible accelera¬ 
tion, o f star forma tion that would have happened anyway 
(IDale et al.ll2007bll . However, in feedback dominated envi¬ 
ronments, back-pressure from the photoevaporative flow off 
the globule will continue to compress the globule even as 
the protostar evolves. From detailed study of the ‘de fiant fin- 
;er,’ another small globule in the Carina nebula JSmith et al.l 
2004bli find that back-pressure from the photoevaporative 
flow is a factor of ~ 5 greater than thermal support in the 
globule. If similar conditions apply to HH 900, then pressure 
from the ionization front may alter the physical conditions of 
the local star forming environment. By increasing the local 
pressure, feedback from massive stars may affect the enve¬ 
lope accretion rate, the amount of molecular and dusty ma¬ 
terial entrained by the jet (as well as the degree that it can 
clear an outflow cavity in the molecular envelope), and the 
rate at which the developing protostar destroys the globule 
from within. Clearly, hydrodynamic models of the evolution 
of protostars in irradiated globules would be interesting. 


5 CONCLUSIONS 

We present new imaging and spectroscopy of HH 900, a 
peculiar protostellar outflow in the Carina nebula. Ha im¬ 
ages reveal an unusually wide-body bipolar outflow while 
new narrowband [Fe 11 ] 1.26 pm and 1.64 pm images of 
HH 900 obtained with WFC3-IR trace a collimated jet. 
[Fe 11 ] emission remains remarkably symmetric on either side 
of the globule, suggesting that the unseen driving source lies 
roughly at the center of the globule. New narrowband H 2 im¬ 
ages from GSAOI reveal extended molecular emission from 
the outflow just outside the globule that disappears at the 
same separation where the [Fe 11 ] emission from the jet grows 
bright (~ 1.5"). The H 2 —> [Fe 11 ] transition zone is seen on 
both sides of the outflow and is roughly symmetric. Both 
Ha and H 2 emission are as wide as the globule when they 
emerge into the H 11 region. Their width and kinematics 
suggest that these lines trace dusty and molecular material 
entrained by the jet that was recently dragged into the H 11 
region. 

Optical and near-IR spectra reveal separate and distinct 
velocity components traced by [Fe 11 ], H 2 , and Ha emis¬ 
sion. Both [Fe 11 ] lines trace the steady jet emission from 
the red- and blueshifted (eastern and western, respectively) 
limbs of the jet. On the other hand, velocities in the Ha and 
H 2 spectra increase further away from the globule, similar 
to the Hubblc-like velocity structure observed in molecular 
outflows. 

Unlike other HH jets with a S/nteer-identihed candidate 
driving source near the jet axis, [Fe 11 ] emission from HH 900 
does not connect the jet to either of the two YSOs near 
the globule. The jet axis defined by [Fe 11 ] emission clearly 
bisects the globule, and dust lanes near the globule surface 
reveal the impact of the jet breaking out of the globule. 
Together, this argues for the existence of a third protostar 
that remains deeply embedded in the globule and is obscured 
by high column density, perhaps by an edge-on circumstellar 
disc. The invisibility of this protostar may be interpreted 
as evidence that the source is extremely young, that the 
globule has been compressed to high densities by radiation 
from nearby massive stars, or both. 

Our data allow us to reject the hypothesis that the YSO 
(PCYC 838) along the western limb of broader bipolar out¬ 
flow drives a separate microjet. The [Fe 11 ] jet structure is 
symmetric and center on the globule (not the star) and there 
is no evidence for multiple velocity components in either Ha 
or [Fe 11 ] spectra. Proper motions measured in Ha images 
taken with HST ~ 9 yr apart are consistent with zero veloc¬ 
ity, calling into question the nature of this unusual feature. 
Low velocities rule out a dynamical origin. Nevertheless, the 
small projected separation between HH 900 and this YSO 
suggests that the YSO may illuminate the nearby bipolar 
outflow, creating an Ha-bright filament. However, the sim¬ 
plest interpretation is that this YSO is simply a chance align¬ 
ment with the HH 900 outflow. 

The bright, collimated [Fe 11 ] jet, together with wider- 
angle Ha and H 2 components point to a powerful jet break¬ 
ing out of a small globule that is sufficiently massive to have 
survived the harsh radiative environment of the Carina neb¬ 
ula. Different morphologies and velocity structures from the 
optical and IR emission lines in HH 900 reveal the protostel¬ 
lar jet itself and offer a rare glimpse of the material it en- 
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trains before it is destroyed in the H II region. This supports 
the interpretation that the irradiated outflows in the Carina 
reflect the same underlying phenomena driving molecular 
outflows in embedded regions, but that they are laid bare 
and illuminated by the harsh radiative environment. 
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